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Abstract 



Topological excitations in Chern-Simons (CS) gauge theories which describe double- 
layer fractional quantum Hall effect are studied. We shall consider generic (m, m, n) 
Halperin state. There are two type of solitons; one is vortex type excitation which has 
essentially the same structure with quasi-hole excitation in the single-layer case. The 
other is nontrivial pseudospin texture which is so-called skyrmion or meron. We shall 
first study qualitative properties of solitons in the original CS gauge theory and give 
results of numerical calculation. Then by using duality transformation, we derive an 
effective theory for topological excitations in the fractional quantum Hall effect. For 
spin texture, that theory is nonrelativistic CP 1 nonlinear a-model with a CS gauge 
interaction and a Hopf term. Finally, we study quantum mechanical system of multi- 
soliton states, retaining only center coordinates of solitons as collective coordinates. 
Existence of inter-layer tunneling drasically changes excitation spectrum. 
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1 Introduction 



Study of 2-dimensional electron system in a strong magnetic field is one of the most 
interesting topics in the condensed matter physics. Especially in the last few years, 
both experimentally and theoretically, existence of topologically stable excitations is 
suggested for quantum Hall state (QHS) with internal degrees of freedom like spin 
singlet QHS and QHS in double-layer systems. Theoretical study started with QHS at 
filling factor v — 1 by taking into account spin degrees of freedom]]]]. It was suggested 
that low-energy charged excitation in that system is not a single-body electron but a 
topological spin texture with finite size and that excitation has similar structure with 
skyrmion in ferromagnet. At present, some experiments support this observation [[J. 

Somewhat detailed study on the double-layer case is given in Ref. by using lowest 
Landau level (LLL) projection in the first-quantization language. Most of studies on 
double-layer QHS given so far are restricted to Halperin's (m, m, m) states. 

In this paper, we shall study topological excitations in generic (m, m, n) states in 
the framework of the Chern-Simons (CS) gauge theory. In the CS gauge theory for 
QH effect, the original fermionic electron is expressed by bosonized electron with odd 
number of CS flux quanta attached^]. QHS is characterized by the coherent conden- 
sation of bosonized electron. In the single layer case, there is a topological excitation, 
which has a close resemblance to the vortex in the superfluid, and that excitation is 
identified with Laughlin's quasi-hole excitation in the fractional quantum Hall (FQH) 
state 0. In the QHS with an internal symmetry, other topological excitations are 
expected which have nontrivial structure in the internal space. 

This paper is organized as follows. In Sec. 2, we shall explain models, which 
describe FQH effect (FQHE) in double layer (DL) electron systems. In general, 
there are two CS gauge fields corresponding to U{1) ® U(l) local symmetry, and 
two bosonized electron fields, which represent electrons in upper and lower layers, 
respectively. Ground state at specific filling factor has coherent condensation of both 
bosonized electron fields. In Sec. 3, it is explained that there are essentially two types 
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of topological solitons in the models. One is vortex which has nontrivial winding 
number for U(l) symmetry corresponding to the simultaneous phase rotation of the 
pseudospin up and down bosonized electron fields. The other is meron which has 
nontrivial topological structure in the SU(2) internal space of pseudospin. Qualita- 
tive study of the topological solitons and numerical solutions to field equations are 
given. In Sec. 4, an effective field theory for topological excitations is derived by using 
duality transformation. We first decompose bosonized electron fields into amplitude 
or density, common U(l) phase factor and CP 1 variable. Lowest Landau level condi- 
tion is obtained explicitly, which determines deviation of charge density in terms of 
topological currents of vortex and meron. Pseudospin CP 1 part of the effective field 
theory is a CP 1 nonlinear a- model interacting with a CS gauge field and with a Hopf 
term. In Sec. 5, quantum mechanics of solitons is discussed. Only center coordinates 
of them are taken into account as collective coordinates. Effect of inter-layer tunnel- 
ing is considered. Single meron cannot exist with finite energy. Spectrum of meron 
pair is studied. Section 6 is devoted for discussion and conclusion. 



2 Models 

Lagrangian of the model is given as follows, 
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E \D°M 2 ~Vm, (2-1) 



2M 

c=T,J. 

£ cs = £ csK) + £ csK) 

= -\^ux{^a+d v at + -^a-d v al), (2.2) 

where tp a (a — 1, 2 or f, J,) is the bosonized electron fields in upper and lower layers, 
respectively, M is the mass of electrons and p and q are parameters. Greek indices 
take 0, 1 and 2, while roman indices take 1 and 2. e^ u \ is the antisymmetric tensor 
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and the covariant derivative is defined as 

Df = dj - iaj =f iaj + ieAj, (2.3) 

where external magnetic field is directed to the z-axis and in the symmetric Coulomb 
gauge Aj = —^€jkXk- V[i^ a ] is interaction term between electrons like the Coulomb 
repulsion and short-range four-body interaction, e.g., 

5i(|^tI 4 + I^I 4 )+^ 2 I^I 2 ' etc - (2-4) 

Chern-Simons constraints are obtained by differentiating the Lagrangian with respect 

tO CLq 

eijdidj = 2q^a 3 ^ 

* = (V>T ^if, (2.5) 

where = eoij- As we shall see, (p + q) must be an odd integer times 7r and (p — q) 
must be an integer times n for the original electrons to be fermionic. 
Partition function in the imaginary-time formalism is given as 

Z = J [D^D^Da] exp{ J dTd 2 x(C^ + £cs)}- ( 2 - 6 ) 
The Lagrangian C = + £q§ is invariant under U(V)(8)U (1) gauge transformation, 

\J/ — > g»01+»02<T3\J/ 



a 



< + d ll 6 1 (2.7) 



a M -> a M + <9 M # 2 . 

Ground state for fractional quantum Hall effect is given by the following field 
configuration, 



eB = eydidjo = 2pp, e^a^ = 0, a+ = eAj, (2.8) 

v =-£ = -, 2.9 
eB p 

where p is the average electron density. It is easily verified that the above static 
configuration is a solution to the field equations if and only if the filling factor v has 
specific values given by (|2.9|) . 

One may think that there is a massless mode, because the configuration ( |2.8| ) 
breakes the U(l) <g> U(l) symmetry. However, these massless modes are eaten by the 
CS gauge fields and the CS gauge fields acquire mass term by the Higgs mechanism. 
This is important, because QHS must be incompressible and there must be no gap less 
modes. 

It is not so difficult to show that the above ground state (|2.8| ) corresponds to the 
following Halperin's (m,m,n) state|6j in the first-quantization language]?]; 

= ri(4 -z]r n(4 -*jr n(4 -4r 

i<j i<j i,j 

eB 



x 



^(-x^l' + Ekl 2 )), (2-10) 



where m = ^(p + q) and n = ~(p — q). For the specific case m = n, one can set 
aj = and only the single CS gauge field remains. 



3 Solitons and numerical calculations 

In this section, we shall study topological solitons in the CS gauge theory ( |2.1|) and 
( p.2| ). We first assume the following specific form of the potential between bosonized 
electrons; 

%] = + 1^3*) 2 + wm, (3.i) 

where W^Vv] represents long-range inter-layer and intra-layer Coulomb repulsions. 
It is argued that the above short-range repulsion is required in order to retain the 



Pauli exclusion principle in the bosonization of fermionic field ||. In the present case, 
field equation to be solved becomes very simple with this term as we shall see. It 
is expected, however, that the existence of this term influences behavior of solutions 
only near the center of them. 

It is easily verified that the following modified Bogomol'nyi decomposition holds, 

a a a 

- e^(a+ - eAj) - e^aj (#<r 3 tf) . (3.2) 

By using the CS constraints (|2.5|) and (|3.1| ), Hamiltonian (in the real-time formalism) 
is given as 

w = 2^E \( D i - iD M + y + WW.], (3.3) 

where u c = eB/M. If we treat the Coulomb interaction as a perturbation, the lowest- 
energy configurations satisfy the following "self-dual" equations, 

(£>J - iD\)^ = 0, {D[ - iDi)^ = 0. (3.4) 

Actually, (|3.4j) is nothing but the lowest Landau level condition, and there are many 
solutions to ( |3.4| ). This degeneracy is solved by the Coulomb interaction, and the 
ground state is the uniform configuration ( j2.8| ). 

We are going to seek topological soliton solutions to the self-dual equation ( |3.4|) . 
In this section, we are interested in single soliton configurations. Therefore, we pa- 
rameterize the fields as 

V>f = V^exp (ui(r) + in 1 9{x) s j 1 

i*l = v^exp (u 2 (r) + in 2 9(xf), (3.5) 

where r = \x\, 6(x) is the azimuthal function 9{x) = arctan and n\ and n 2 are 
integers which label topological solitons. By substituting ( |3.5|) into (|3.4| ), we have 

a t + a k - eA k = n x d k 9 - ChAm, 

a\ - a k - eA k = n 2 d k 9 - e ki diu 2 . (3.6) 
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By differentiating Eq.( |3.6| ) and using the CS constraint ( |2.5| ), we have the following 
Toda-type equations for Ui{r) {i = 1,2), 

¥ + AdA +ni T + 2 -^^ e =0 ' (3 ' 7) 
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r- l' 

X = VWj -'^Wo' K = P 



(3.8) 



p + q p — q 
yp - q p + q 

where Iq is the magnetic length l = 1/y/eB. It is obvious that u\{r) and w 2 (r) couple 
with each other nontrivially unless p = q. 

It is useful to rewrite ( |3.7| ) by introducing t>j(A) = Ui(X) + n^biA as follows; 

(^ + ^- + 2 - 2 W' A ^ - < 3 - 9 » 

where we have used the identity 
In terms of fj(A), 

tpi = \/pexp (vi(X) — riilnX + ini6(x)^j, 

i>l — y/pexp (v 2 (X) - n 2 \n\ + in 2 9(x)y (3.11) 



From (|3.11 ), it is obvious that the integers rii must be nonpositive, rij = 0, —1, —2, 



for single- valuedness and regularity of solution at the center of soliton, and t>j(A)'s are 
regular functions. On the other hand, boundary condition at the infinity is given as 

lim e Ul{x) = -^=, (3.12) 

which guarantees that at sufficiently far from the soliton, amplitude of solution is 
equal to that of the ground state. 

Before going into detailed study of the numerical calculation, let us examine qual- 
itative properties of solitons. There are essentially two types of topological solitons. 
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The first one has a nontrivial winding number for £7(1) symmetry of simultaneous 
phase rotation of the both up and down bosonized electron fields. This soliton corre- 



sponds to the (ni, n 2 ) = (—1, —1) in (|3.5|) , and has trivial structure for the pseudospin 
internal space. We shall call this soliton vortex. As we shall show in the numerical 
calculation, behavior of this soliton is very similar to the vortex in the single-layer 
system. The other ones are pseudospin texture and correspond to the (—1,0) and 
(0, —1) in (|3.5| ). Functions wi(A) and «2(A) for (—1,0) behave as 

«i (A) ~ In A, «2(A) ~ const., for A — > 0, 

Ul (X) = u 2 {X) \nV2, for A ^00, (3.13) 

and similarly for the (0, —1). It is useful to introduce "pseudospin operator" m(x) as 
usual byQ 

fh = ^a^. (3.14) 



From (|3.13| ), it is obvious for the (—1, 0) meron 
m(x) ~ — z, as x — > 0, 

fh(x) ~ cos#(x) • x + sin6'(x) • y, as |x| — > 00, (3.15) 

where x, etc., are unit vectors in the x, y and z directions. Picture of the above 
configuration of rh(x) is given in Fig.l. Similar picture is obtained for (0, —1) soliton; 

fh(x) ~ z, as i^0, 

rh(x) ~ cos6'(x) • x — sin#(a;) • y, as |x| — > 00. (3.16) 

It is obvious that they have fractional winding number in the 0(3)-pseudospin space, 
and we call them meron. A pair of (—1, 0) and (0, —1) merons form a skyrmion with 
unit topological number.^ 



1 rh(x) is not invariant under transformation ( |2.7| ). But under regular gauge transformation, 
topological structure of rh{x) is invariant. Gauge-invariant definition of the topological charge will 
be given in Sec. 4. 

2 More precisely, we should introduce pseudospin operator with unit length m = rn/\rh\. Then, 
topological charge density is given as q{x) = -g^e^m ■ [dim x djih}. 



Let us introduce a couple of quantities; soliton charge Q is denned as 



Q = J d 2 x(^ty - ^o*( 



— fdx^af-aloj 



cs 



2p 



(3.17) 



where \l/o etc., is the quantities of the ground state given by ( |2.8|) , and $ cs is the total 
CS flux of soliton. The electric charge Qe is Qe = — eQ- Similarly we define 



$ cs = j <7,r' a 

Q 



cs 



2q 



Total angular momentum is given by 



1 

27 



J = /few, 

^T-DjVt - BJ*t • </>t} + {^Bjj/'j - BjiAj • 4>i} 

It is useful to decompose the angular momentum as J = U + V, 
1/ = -eij J d 2 x Xi(af^^ - aj^a 3 ^ 



The CS constraint Q2.5Q is solved in the Coulomb gauge as follows, 



a\(x) 



a- [x) 



P 

7T 



F — 2/1 
(x - y) fc 



F - 2/1 



*(r 3 ^(y). 



From (ETZ 



(3.18) 



(3.19) 
(3.20) 



(3.21) 



(3.22) 



(3.23) 
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On the other hand for U, from ( |3.5| ) and ( |3.6| ) 

U = Y,n a [d 2 x \ip a \ 2 . 
Total angular momentum of soliton, J(soliton), is naturally defined as 



(3.24) 



J(soliton) = J[ip a ) - J[ip a: o = \J p/2] 

= e^/ <?*m 2 -^[{J d 2 x m} 2 - {Jd\ wof 

~[{ J d 2 x ^o^} 2 -{J ' d 2 x ^s^o} 2 " • 
Similarly, energy of soliton is given as 

E = Jd 2 x{^Y.\( D i-* D 2)^\ 2 + Y^ + W ^° 
- J d 2 x{^ ^ + W[ij afi }} 
'd 2 x{w[^]-W[^ }}, 



(3.25) 



(3.26) 



where we have used the self-dual equation (3.4) and the fact that the total number of 
electrons is conserved. For general configuration iJj a (r)e m,T9 of the winding numbers 
(rii,n 2 ) and the CS gauge fields of the following form 



at = dh.9 ■ a ± ('r), ar = at ± a h = dk.6 • a u (r 



h u k 



Hi 



(3.27) 



the energy is explicitly given as 



E 



2vr / — 

r 



dr r 1 



2M 



E 



# CT \2 



dr 



eB 



dr 



+ ( a ° - e Ji r 2 - n ff ) 2 (i) 2 } + ^r 2 (M - * * ) + r 2 (w[^] - W$ aA 



2 ' yr ' J 2 

Therefore, for energy to be finite, the CS gauge fields must satisfy^ 

eB 1 
lim a + (r) = —r 2 + ~(ni + n 2 ), 

r— >oo 2 2 

lim a~(r) = -(rii — n^)- 

r^oo 2 



(3.28) 



3 For the (m,m,m) state, there is only single CS gauge field a+ and = 0. Then, merons have 
logarithmically divergent energy H. 



10 



From (ggg ), 

Q = ^(^). o = / - w.) = 

The above J(soliton) should be interpreted as spin of solitons, S. On the other hand, 
statistical parameter of soliton a (soliton) is given as follows by the argument of the 
Aharonov-Bohm effect, 

a(soliton) = a + + ct_, 

«+ = — y~ = ' 

a_ = -9** = - q Q\ (3.30) 

Therefore, the spin-statistics relation, S = — holds generally for solitons in the 
present system. Actually, it is obvious from the above consideration that there are 
two CS gauge fields, a+ and a~ , in the present model and they contribute to both spin 
and statistics additively. This result is rederived in later discussion in Sec. 4 which is 
based on an effective theory of solitons. 

Let us turn to the numerical solutions to the field equations (|3.9|). For example 



for the (—1,-1) vortex, asymptotic solution near the center of soliton is obtained as 

Wi (A)~A-y, (3.31) 

where /3j's are parameters to be determined by the boundary condition at spatial 
infinity. Similar asymptotic solutions are obtained for general (n 1 ,n 2 ). Using these 
asymptotic solutions, we obtain numerical solutions. 

In the DL FQHE, there are more than one QHS's at the same filling factor. 
Actually as we showed in Sec. 2, the present CS gauge theory describes QHS at v — 
tt/p, which is independent of the parameter q. The value of q is determined by the inter 
and intra-layer Coulomb interactions W^], which depend on inter-layer distance d. 
For each d, some specific Halperin state is a good variational wave function. 
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As an example, we shall consider the case v = |, and at this filling fraction, 
Halperin's states (3,3,2), (5,5,0) and (1,1,4) are possible candidates for the QHS. 
From (|2.10| ), it is obvious that the state (5,5,0) corresponds to such a case that the 
inter-layer distance d is sufficiently large and the inter-layer Coulomb interaction is 
negligibly small. On the other hand for the (1,1,4) state, the inter-layer Coulomb 
interaction dominates over the intra-layer interaction. Though this situation does 
not seem to be realized in samples in experiment so far, some interesting features of 
solitons appear in this QHS as we shall see. The state (3, 3, 2) is in between. 

Vortex solution {n\ = —l,n2 = —1) in the (3,3,2) state is given in Fig. 2. The 
picture shows that size of vortex is about four times of the magnetic length l = 
l/\/eB, and its form is quite similar to that in the single-layer system ||. Meron 
solutions (ni = — 1, n-i = 0) etc. are also given in Fig.2. It should be remarked that 
the density in the lower layer increases at the center of the meron. It is a result of the 
inter-layer Coulomb interaction. Size of the meron is about 4.5\/2Zo, and it is larger 
than the vortex. Creation energies of vortex and merons are calculation from these 
solutions and ( |3.26| ). We also show the numerical solutions of (—2,0) and (—2,-1) 
solitons in Fig.2. 

Similarly, meron solution for (5, 5, 0) state is shown in Fig. 3. As there is no inter- 
layer interaction, the density in the lower layer is not influenced by the deviation of 
density in the lower layer. Finally for the (1, 1,4) state, meron solution in Fig.4 has 
a long undulating tail because of the dominant inter-layer interaction. 

From the above meron solutions, we can calculate normalized pseudospin vector 
rh(x) = m{x) /\m(x)\. Merons have a fractional topological charge in the 0(3) space. 
In the following section, we shall derive an effective field theory of solitons, and show 
which field theory describes rh(x) and how the topological charge is related with 
electric charge of solitons. As we see, the lowest Landau level projection plays a very 
important role there. 

Finally, we shall give wave functions corresponding to the topological solitons in 
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the first-quantization language. For soliton of (rii, n 2 ) type which is located at (w, w), 



MM**) = ri(4 -4r ri(4 -4r ri(4 -4r 



xIl(4-^)- m II(4 



— u> 



X 



ex P (-^(Ei4i 2 + El4i 2 )). 



(3.32) 



4 Effective action of vortices and merons 

In this section, we shall derive an effective field theory of topological solitons con- 
sidered in Sec. 3. At specific filling factor ( |2.9|) , these solitons are excitations of low 
energies. Lowest Landau level condition is expressed in terms of topological current 
of these solitons, as we shall see in later discussion. 

We first prameterize the bosonized electron fields as follows ^] , 



e fe <p v e U{1), z = { Zl z 2 y e CP 1 , 



(4.1) 



where e ld is a regular part and <fr v is a topologically nontrivial part, which represents 
vortex degrees of freedom. CP 1 variable z a represents pseudospin degrees of freedom 
and satisfies CP 1 condition \ z a\ 2 — L and ip a is invariant under z a — > e ta z a and 



Substituting (Q into Q), 



where 



-Jo (pdo4> + (zdoz) — iog — ia Q {za^z) 

" 1 + zdjZ — i(aj~ — eAj) — iaj(za 3 z) 

dj(j) ■ <p + djZ ■ z + z(a+ — eAj) + iaj (za^z) + C z 
Jo 



2M 



2M 



Djz ■ DjZ + (zDjz) 



D^z = (dp - ia^a 3 )z, 
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(4.2) 



(4.3) 



and we have neglected derivative terms of Jo, for Jo is not a dynamical variable^ and 
its value is determined by LLL condition as we shall see. We have also neglected the 
interaction term V^[Vvl> which will be discussed later on. 

We shall perform duality transformation [|TT| . First, we use the following identity; 



In the present case, K = and 



+ z ■ —z 

% ' i 



(aj - eAj) - a j (zcr 3 z), 



(4.4) 



(4.5) 



and the Lagrangian C is expressed as 



C 



-Jn 



%cj) + (zd z) - ia,Q - ia (za 3 z) 



2K J - J ~ U > 



)—<p + z ■ —z 

i i 



+iJj(aj - eAj) + iJjcij (za 3 z) + C z + Cq$. 



(4.6) 



It is obvious that integration over 9 gives the divergentless condition 



<V M = o. 



(4.7) 



Solution to 



is 



Jn — —ttwxdybx. 

Z7T 



Substituting into fl£E|) , 



AH 
1 1 



4>vd 4>v + (zd z) - ia^ - ia (20-32) 



2if (2tt) 2 



(<9 &j - ^6 ) 2 - 7^£ij(d bi - dib 0j 



Z 9 
I 



-& v + z ■ —z 



+— e ij (d b i - dib )(aj - eAj + aj(za 3 z)) + C z + £ cs . 



(4.8) 



(4.9) 



Next we integrate over 60 with the Coulomb gauge djbj = 0. The form of the 
integrand is the following, 

1 1 



2 A -(2,)i (9 ' i '» )2 + !i '»*' 



-—end, 



2tt 



5, 



'v^r^v + z ■ —z - (at - eA- 



a,- Uo- 3 zl 



(4.10) 



4 It is obvious that the term Jq^ 2 <9o Jq^ 2 is total derivative 
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Integration gives 

J [Db ]ef d3 ^- c ^ 2+ib ^ = e"*/*^" 1 *, (4.11) 



where C = and 



J $A _1 $ = J dr J d 2 xd 2 y $(x) In \x - y\$(y), (4.12) 
We furthermore integrate over the CS fields Oq which gives the CS constraint, 



^fiat = 2p\^\ 2 = 2pJ = -eijdibj, (4.14) 

7T 



and $ is expressed as 



1 v -, v % v + z.Rl z -(P b eAj)]. (4.15) 

2 



Let us consider quantization of the field bj. Kinetic term is given as (dobj) 
and a system is simply a harmonic oscillator with the following form of the action 
L = x 2 + (x — a) 2 . We shall consider static configuration of solitons. Then, Largangian 
of bj field, Cb, is given by 

C b = -C(d b J ) 2 -^J d 2 y$(x,T)A- 1 (x-y)$(y,T). (4.16) 

It is convenient to take B(x) = tijd^bj as dynamical variable for quantization in stead 
of bj(x) itself which is subject to the Coulomb gauge condition. The original field bj 
is expressed in terms of B as 

h{x) = e l3 d] J d 2 yA'\x - y)B{y). (4.17) 

We shall introduce Fourier transformed variables as 



Ms) = l^ qX ~ h Ml 

B{x) = J** e *°B(q), etc. (4.18) 
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In terms of Fourier transformed fields, 

d 2 q 1 



d 2 x (dobjY 



(2vr) 2 g 



2 d B(q)d B(-q) : 



(4.19) 



Then, returning to the real-time formalism as t = —ir, canonical conjugate variable 
of B(q) is 



nB(9) = C (2W^ 9) - 



Hamiltonian is given as 



K 



K ~, s 1 ~ 



and 



Hb = / d'q -q i (2^)'n B (q)n B (-q) + -*(,)-*(-,) 



(4.20) 



(4.21) 



(4.22) 



From ( 4.21 ) and (^4.22 ), it is obvious that the ground state is the configuration which 
satisfies $ = and excitation energies are given by Nuj c , where N is the number of 
excitations and u c = eB/M. 

Ground state condition $ = is nothing but the LLL condition, which is given as 
follows 0, 

Jo + Jo - + |^A, = 0, (4.23) 

where topological densities and currents are defined as 

Ox 



Jl = ^^A^(^y^), 

4 = —e^ x d v (z^z). 



2vr 
1 

2tt'' v z 

There is another U(l) (g) Z7 (1) gauge-invariant topological current, which is defined as 



(4.24) 
(4.25) 



J u = —^nu\Ou\ZU^—Z 
271 x I 



(4.26) 



where D\ = d\ — i{a\ — eA\). This topological charge distinguishes between (—1, 0) 
and (0, —1) merons. 
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In the small-size or the long-distance limit of the topological solitons, the above 
currents are given as 

a 
a 

Jo(x) = J2q S a S(x-x s ' a ), 

a 

Jf{x) = "£q S a xf a 5(x-x s > a ), ( 4 - 27 ) 

a 

where index a distinguishes solitons and q v a and gf are topological charges of solitons, 
i.e., 

q v = —1 for (—1, —1) vortex, 

q s = — - for (—1, 0), (0, —1) merons, 

q s = +^ for (1, 0), (0, 1) merons, (4.28) 



and x v ' a and x s,a are their coordinates. From ( |3.11|) , merons (1,0) etc. do not exist 



as a regular solution to the self-dual equation. This situation is similar to the single- 
layer system, in which quasi-particle excitation does not exist Kular solution. 
However, these solitons must exist as low-energy excitations and their behavior away 
from the center is described by ( |3 . 1 1| ) . 



The LLL condition (f4.23|) determines the density of electrons in terms of the 
topological currents. For example from ( |4.14| ), the vortex (—1, —1) has charge ve and 



the merons (—1,0) and (0, —1) have charge ^. 



Here, a few remarks about the LLL condition are in order. The condition ( [4.23 



is 



correct for static configurations. For moving soliton, it is not so difficult to show that 
there appears additional term which is proportional to e^S^JJ + Jj). Existence of 
this term, however, does not change the discussion on statistics of solitons, etc, given 
below. For example, 

/ d 2 x eijdiJj = / Jjdxj, (4.29) 

J So J dSo 
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where So is a domain and OSq is its boundary. For sufficiently large domain whose 
boundary is far away from solitons, the above integral is vanishingly small and there- 
fore charge of solitons is determined solely by the topological density. 

There is obviously an ambiguity in the decomposition (O), i.e., the bosonized 
electron fields are invariant under the following transformation, 



-up 



', Z a > 6 ^Zfj, 



(4.30) 



where the phase <p can have a nontrivial winding number. This singular transforma- 
tion is allowed at the center of solitons with topological numbers rii ^ and n 2 ^ 0, 
as the amplitude Jo vanishes there. From the definition of the topological charges 
( Oj ) and flOg) , 



Z7l 



(4.31) 



Therefore, the total topological charge is invariant under ( |4.30| ). 



With the LLL condition ( [4.23D , the Lagrangian is given as 



£ = -—tindib 



2tt 



— i i r — ■ d ■ 

4>vd <f)v + {zd z) - icio ■ {za 3 z) - —eijOoh 4> v —<f) v + z ■ —z 

J 2tt 1 i i 



% % 
+ 7r e ijdobi ■ (a| - eAj) + —e^dobi a~ ■ {za 3 z) + C z + £ cs (aC). (4.32) 

After some calculation, it is shown 



lst+2nd+4th terms of (4.32) = -ify ■ (J? + J? ). 



(4.33) 



On the other hand, the third term is evaluated as follows; 



—eijdok ■ (a+ - eAj] 
In J 



—eijbi ■ <9 (a+ - eAj) 
In J 



Let us denote J M = J v + . In the Coulomb gauge, 



(4.34) 



(at - eAj){x) = 2vr / <Pye jk d£A-\x - y) J (y). 



(4.35) 



By the conservation of the topological currents 

d^Jf, = 0, 



(4.36) 



6ijd (aj - eAj)(x) 



-2tt J d 2 yd*A- 1 (x-y)d j (y) 

2vr J fPyd&A-^x - y) ■ J k (y). (4.37) 



On the other hand, 



7 / \ 2tt 2 , , 
bi(x) - 1 



P 



71 



d^ye^A 1 (x-y)- J (y) + -eA^x 



V 



(4.38) 



Then from ( [OTp and QO^) 



y d 2 xbi(x)eijdo(a^ — eAj)(x) 



2tt2 

p 



d 2 xd 2 yd 2 z 



d?d%A-\x - z) ■ J k (z)e^A~\x - 



7T 



xJ (y) + - / d^ze^faOfiffljfA-Va; - z) • J fc (z) 
J p 



= 0, 

where we have used the Coulomb gauge condition. 

Therefore, the effective Lagrangian Le has the following final form; 

C E = C z + £(js(a~) - ibi ■ Ji. 



(4- 



(4.40) 



With the LLL condition (|4.23 ), the term (|4.33|) is nothing but the Hopf term which 



determines statistics of solitons|T2[. For vortex with (—1,-1), JJ t has charge —1. 
It is not difficult to show that statistical parameter of that excitation coming from 
the above Hopf term is iris, where the filling factor v is given by v = ^ = ^. On 
the other hand for merons (—1,0) and (0,-1), J~ has charge —1/2, and then the 
statistical parameter from the Hopf term is 7rz/(|) 2 . The above results are consistent 
with the discussion in terms of the Aharonov-Bohm effect in Sec. 3, i.e., the Hopf term 
in ( |4.40|) gives the statistical parameter a + in ( |3.30| ). 

The above Hopf term also gives spin to solitons[l2~], which is consistent with the 
spin-statistics relation. As we showed in Sec. 3, merons in the present system have 
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the extra terms both in their spin and statistical parameter, besides those coming 
from the Hopf term. These terms come from the coupling with the CS field a~ , as we 
can see from (|4.40|) and the discussion in Sec. 3. Actually recently, CS gauge theory, 
which is a relativistic version of (|4.40 ) without Hopf term, is studied [13]. There, it is 
shown that solitons similar to those in the present system have a fractional spin and 
also fractional statistics. We also studied topological solitons in the effective theory 
Q4.40Q very recently|14| and found that the CS gauge field a~ gives spin and statistical 
parameter to merons which are consistent with the results in Sec. 3. 

In the following section, we shall consider the quantization of multi-soliton con- 
figurations retaining only center coordinates of solitons. In that discussion, it will 
become clear that the above Hopf term surely determines statistics of solitons. 



5 Quantum mechanics of topological excitations 

In the previous section, we have obtained the effective Lagrangian of topological 
solitons Q4.40D by the duality transformation. Also as we showed by the numerical 
calculation in Sec. 3, topological solitons have a finite size of order of the magnetic 
length. Actually form of the meron is determined by the terms C z + £(^g(a~) in C E , 
and long-wave-length properties of topological excitations are determined by the term 
6 j • Jj. Especially, by using the topological currents in the form of the small-size limit 
( 4.271 ), one can easily obtain an effective Largangian of multi-soliton system from 



hi ■ J{, which retaines only coordinates of solitons as collective coordinates. 

The LLL condition is solved as ( j4.38|) . It is not so difficult to calculate the long- 



distance limit of the term hi ■ J t , 



J d 2 xbi(x) ■ Ji{x) = lipEij ( C^T^T + Yl VaXi 



/ v v * !l i V s v ■S,a\' L x J J 



r .,, _ x v,/3)2 ~ Z-^ HaH^ij-^i ( x S,a _ x v,/3)2 
I a. 8 V / 
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i v S ■ v 




4 + Qa1p e ij X i 




(5.1) 



It is obvious that the terms proportional to v in ( |5.1|) determine the statistics of 



where 9{x) = arctan (jj^j- In the path-integral formalism, this statistical term gives 
nontrivial contribution when world lines of two particles entwine each other. As 
we briefly mentioned in Sec. 2 and Sec. 3, similar terms appear from the coupling of 
solitons with the CS gauge field a~. 

Anyonic properties of topological excitations plays an important role, for exam- 
ple, for hierarchical structure of the FQHE. Anyonic vortices and merons can be 
expressed in terms of boson fields with CS flux attached. At specific filling factors, 
these bosonized anyons can Bose condense, as the bosonized electrons. These states 
correspond to daughter states of the FQHE, and it is interesting to compare theoret- 
ical prediction of filling fractions of the FQHS with experimental observations ||15|| . 

In this section, we shall study quantum system of quasi-excitations in the FHQS. 
The first term of ( |5.1|) determines canonical commutation relations (CCR) of soliton 
coordinates. Generally, the quantum system of anyons, whose Lagrangian is given by 
QS.ip , is very difficult to solve. Therefore, we shall mostly concentrate on two-particle 
systems in the rest of discussion. In specific but physically interesting cases of two- 
particle systems, it can be shown that the statistical terms can be ignored in the first 
approximation, as we discuss later on. On the other hand, we shall take into account 
interactions between solitons. 

Returnning to the real-time formalism, the Lagrangian of the multi-soliton system 
5 It should be remarked that in the present system, the size of solitons or the magnetic length 
itself plays a role of short-distance cutoff. 



solitons, because they are rewritten asf] 




(5.2) 
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is given as, 

£ soliton = -npeij{±^Qa x i x j +2^9^ 



S,a S,at 



a lu,lu'=v,S a,/3 

where D(q£) = D v and D(q%) = An are energies to create the vortex and the meron, 
respectively (see Q3.26p ), and W(x UJ ' a — x^ is the Coulomb interactions between 



solitons. From (|5.3| ), the CCR 's among soliton coordinates are given asp6 



\x v,a x"' 13 ] = -i 6ij 8 r 



[xr,xf P } = 0, 



Here we have neglected the statistical terms. After obtaining quantum mechani- 
cal states of two-particle systems, we shall give plausible argument supporting this 
approximation. From (|5.4j ), conjugate momentum of This obviously cor- 

responds to the cyclotron motion in the classical picture. From this fact, dispersion 
relation of a pair of solitons is, in some cases, obtained quite easily. 

Actually concerning vortices, the Lagrangian ( |5.3| ) has a similar structure with 
that of vortices in the single-layer system. It is rather straightforward to show that 
a vortex-antivortex pair, which are separated with each other more than the magetic 



length, form a magneto-roton excitation. Its dispersion relation is given asfTT 



3 2 

E(p t ) = 2D V - (5.5) 

where pt is the eigenvalue of the total momentum. This result is obtained by ignoring 
the statistical terms, but is in qualitative agreement with the results of the exact 



diagonaliz at ion [[171 



Similar relation is obtained for a pair of merons. However, there is a very im- 
portant effect which we have to take into account in the DL electron systems, i.e., 



22 



inter-layer tunneling. Hamiltonian of the inter-layer tunneling is given as, 



= -AsAsfetf), (5.6) 



where Ag^g is the energy splitting between symmetric and anti-symmetric states 
with respect to the pseudospin index.f] 
In terms of the parametrization (|4.1|) , 



H T = -Ag AS J (>(Xi2), (5.7) 

and therefore it is obvious that T~Ct gives an linear attractive confining force between 
(— 1, 0) — (0, —1) pair, (—1, 0) — (1, 0) pair, etc., see Fig. 5. 

Let us first consider a pair of (—1,0) and (1,0) merons, whose coordinates are 
(xi,yi) and (#2,2/2), respectively. From (|5.4j) , the nonvanishing CCR's are 

[x 1 ,y 1 ] = iR, [x 2 ,y 2 } = -iR, R = ~j-- (5.8) 

Introducing the center of mass and relative coordinates as X = \{x\ + x 2 ), Y = 
Myx + 2/2), x = x\ — x 2 and y — yx — 2/2, the CCR' are 

[X, y] = iR, [Y, x) = -%R. (5.9) 



From ( |5.9D , it is obvious that Y = Rp x and y = Rpx, where p x and px are relative 
and total momentum operators, respectively. On the other hand, potential energy of 
this system is give as 

Vl(r 12 ) = -(^) 2 — + T in2 , (5.10) 
v 2 ' r 12 



6 Strictly speaking, the tunneling term (5.6) should be written in terms of the original fermionic 
electron operators. As a result, an extra factor like exp(2z<72 f x ajdxj) appears. This factor is 
important to preserve the gauge symmetry. However, results of the later discussion are not affected 
by this factor. 
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where Tyi = |n — Ti\ and tension T\ oc Ag^g- If the total momentum of the pair is 
directed to the x-axis with magnitude px = Px, the energy is given as 

EM = 2D m - ( y ) w - + ^ Px . (5.11) 

The above dispersion is correct for r 12 > Iq an d Px > lo/R- Therefore, (—1, 0) — (1, 0) 
meron pair also has magneto-roton dispersion. However, in contrast with the vortex- 
antivortex pair, its energy increase linearly with total momentum. 

Let us turn to the (—1, 0) — (0, —1) meron pair. Nontrivial CCR's are 

[x 1 ,y 1 ] = iR, [x 2 ,y2]=iR- (5.12) 

From (1513) , 



[X,Y] = Z -R, [x,y]=2iR. (5.13) 

In terms of the complex coordinates £ = x + iy, [£, !;} = —4R, or by introducing radial 
and angle coordinates as £ = re itp , the CCR is [r 2 , <p] = —4iR. Energy of this system 
is given as 

E(r) = 2D m + (y) 2 i + T 2 r, (5.14) 

where again T 2 oc Ag^g. It is obvious that merons are rotating around the center 
of mass with a fixed radius. The minimum of the energy is located at r m = y^/t^- 
Furthermore at long distance, meron pairs (—1, 0) — (0, —1) and (1, 0) — (0, 1) behave 
as vortices (—1,-1) and (1,1), respectively. 

There are classical pictures corresponding to the above eigenstates of (—1,0) — 
(1, 0) and (—1, 0) — (0, —1) meron pairs. They are given in Fig. 6. Merons are in the 
external magnetic field and also have electric charge propotional to their topological 
charge. Therefore, the Lorentz force works for them, as shown by the first and second 



terms in (|5.3| ). This Lorentz force for merons is balanced by the Coulomb and the 
confining linear potentials between merons in a pair. Therefore, (—1, 0) — (1, 0) meron 
pair moves with a constant velocity and a constant relative position. On the other 
hand in (—1, 0) — (0, —1) pair, each meron rotates around the center of mass coordinate 
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with a constant cyclotron velocity. These classical pictures support the approximation 
which ignores the statistical terms. For (—1, 0) — (1, 0) meron pair, the angle 8(fi — r-j) 
is a constant of motion in the classical limit and then 4i9{f\ — r*2) = 0. In the 
path-intergal formalism, paths close to the classical path dominate. Furthermore, 
the statistical term is invariant under continuous deformations of path. Therefore 
for dominant paths, the statistical term does not contribute to the path integral in 
the present case. We think that this is the reasion why the dispersion relation of 
vortex pair in the single-layer case given by (|5.5|) is, at least, qualitatively correct [ITT . 



Similarly for (—1, 0) — (0, —1) meron pair, 9{ri—r 2 ) oc t in the classical picture, and the 
statistical term gives only irrelevant effects. Therefore quantum states of meron pairs 
which we obtained in the present discussion are consistent with the approximation 
which neglects the statistical terms. We hope that results of exact diagonalization of 
the double-layer systems support out results. 

Here however it should be remarked that in the derivation of the Lagrangian of 
soliton ( |5.3|) , rotation of solitons etc. have not been taken into account, and therefore 
additional terms like (x^' a ) 2 , x^' a x^' a , etc. must appear as higher-derivative terms 
in the effective theory of solitons, if we consider effects of them. For exmple, the 
additional term in the LLL condition eijdiJj induces terms like (x^ ,a ) 2 . This situation 
always happens in deriving an effective theory of topological objects with finite size, 
and most of cases low-energy behavior of these objects is well-described by the low- 
derivative model like (|5.3| ), though many ideas have been proposed for modification of 
that. Especially for constant motion of (—1,0) — (1,0) meron pair, higher- derivative 
terms can be legitimately ignored. 

6 Conclusion 

In this paper, we have discussed topological excitations in the CS gauge theory with 
complex boson fields, which describes DL FQHE. There are two types of topological 
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excitations in the LLL. The one is the vortices which have a nontrivial U(l) winding 
number in the common phase factor of both pseudospin up and down electron fields. 
They are quite similar to the vortex of quasi-hole in the single-layer system. The 
others are the merons which are pseudospin textures and have a nontrivial SU(2) 
winding number. 

We first give qualitative arguments about these excitations and then numerical 
solutions. We derived an effective low-energy theory for them by using duality trans- 
formation and the LLL condition. It is shown that pseudospin textures are described 
by the CP 1 nonlinear a-model with a CS gauge interaction and a Hopf term. The 
Hopf term determines statistical parameters of vortices and merons, and the result 
is consistent with that given by the argument through the Aharonov-Bohm effect. 
We also consider quantum system of multi-soliton configuration by retaining only 
center coordinates of solitons. Lagrangian of this quantum system is obtained from 
the effective low-energy Lagrangian of solitons. We took into account the effects of 
inter-layer tunneling, and determined energy spectrum of meron pairs. It is quite 
interesting to observe these excitations by experiment. Detailed study of dynamical 
properties of topological excitations is under study and results will be reported in 
future publications. 
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Figure Captions 

Fig.l: Configurations of pseudospin vector of merons. Obviously, topological 
charges of them are fractional. 

Fig. 2: Numerical solutions of vortex, merons, etc., in (3,3,2) state of the DL 
FQHE. In both upper and lower layers, charge densities have nontrivial behavior. 
This is a result of inter-layer Coulomb interaction. 

Fig. 3: Meron solution in the (5, 5, 0) state. There are no correlations between 
upper and lower layers. 

Fig. 4: Vortex and merons in the (1,1,4) state. Merons have a long undulating 
tail because of the dominant inter-layer Coulomb interaction. 

Fig. 5: Typical configuration of the (—1, 0) — (0, —1) meron pair in the presence of 
the inter-layer tunneling. In the region between two merons, direction of pseudospins 
is reversed, and there is a confining linear potential between merons. 

Fig. 6: Corresponding classical pictures for systems of meron pairs. Lorentz force 
is balanced by the Coulomb and the linear interactions. 
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(a) (-l,0)meron 




(b) (0,-l)meron 



Figure 1: Configurations of pseudospin vector of merons. Obviously, topological 
charges of them are fractional. 
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Figure 2: Numerical solutions of vortex, merons, etc., in(3,3,2) state of the 
DL FQHE. In both upper and lower layers, charge densities have nontrivial 
behavior. This is a result of inter-layer Coulomb interaction. 
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Figure 3: Meron solution in the (5,5,0) state. There are no correlations between 
upper and lower layers. 
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(a) (-l,-l)vortex (b) (-l,0)meron 

Figure 4: Vortex and merons in the (1,1,4) state. Merons have a long undulating 
tail because of the dominant inter-layer Coulomb interaction. 
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Figure 5: Typical configuration of the (-l,0)-(0,-l) meron pair in the presence 
of the inter-layer tunneling. In the region between two merons, direction of 
pseudospins is reversed, and there is a confining linear potential between merons. 
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Figure 6: Corresponding classical pictures of the systems of merom pairs. 
Lorentz force is balanced by the Coulomb and the linear interaction. 
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